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Abstract: 

Using percolation theory, we derive a conceptual definition of deconfinement in terms 
of cluster formation. The result is readily applicable to infinite volume equilibrium matter 
as well as to finite size pre-equilibrium systems in nuclear collisions. 

1. Introduction 

Quantum chromodynamics predicts that strongly interacting matter will undergo a 
transition from a low density medium consisting of colour-neutral hadrons to a plasma 
of deconfined coloured quarks and gluons at high densities. This transition has been 
confirmed and studied in detail in finite temperature lattice QCD at zero baryon number 
density; nevertheless, even there it can so far be formulated concisely only in the limits 
of infinite and of vanishing quark mass. 

Infinite quark mass leads to pure colour SU(N) gauge theory. Here the Polyakov loop 
L(T) provides an order parameter suitable to define the confinement-deconfinement phase 
transition, driven by the underlying center C SU (N) symmetry of the Lagrangian. 
Confined states share this symmetry, but it is spontaneously broken in the deconfined 
phase. For QCD with quarks of mass m q < oo, the symmetry is broken explicitly in the 
Lagrangian. 

In the limit of vanishing quark mass, the Lagrangian becomes invariant under chiral 
transformations. Now the quark condensate (ipfyiT) constitutes an order parameter, 
distinguishing the hadronic phase of spontaneously broken chiral symmetry from the 
quark-gluon plasma, in which this symmetry is restored. In other words, in the confined 
phase the massless quarks 'dress' themselves to become massive constituent quarks, while 
in the deconfined phase the massless current quarks of the Lagrangian are recovered. Here 
as well, the introduction of a finite quark mass m q explictly breaks the relevant symmetry 
of the Lagrangian. 

The real world lies between these two limits, and there the nature of the transition 
is less clear (Fig. 1). The aim of the present paper is to study this region of broken 

x Talk given at the International Symposium QCD at Finite Baryon Density, April 27 - 30, 1998, 
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1 



symmetries and to look for a more general description of deconfinement which would be 
applicable there as well. Conceptually, the transition from hadronic matter to quark- 
gluon plasma seems rather transparent, no matter what the quark mass is. Once the 
density of constituents becomes so high that several hadrons have considerable overlap, 
it no longer makes sense to partition the underlying quarks into colour-neutral bound 
states. Instead, there appear clusters much larger than hadrons, within which colour is 
not confined. This suggests that deconfinement is related to cluster formation, and since 
that is the central topic of percolation theory, possible connections between percolation 
and deconfinement were discussed already quite some time ago on a rather qualitative 
level 0) 01- in the meantime, however, the interrelation of geometric cluster percolation 
and critical behaviour of thermal systems has become much better understood ||, and 
we want to apply this understanding to clarify the nature of deconfinement. 
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Figure 1: Critical behaviour in finite temperature QCD. 

In section 2, we shall recall some aspects concerning the relation between deconfine- 
ment and chiral symmetry restoration. The main results of percolation theory for lattice 
systems will be summarized in section 3 and extended to continuum percolation in section 
4. This will then provide the basis for the study of percolation in finite temperature and 
finite density QCD in section 5. Finally, in section 6, we shall consider applications to 
deconfinement in nuclear collisions. 

2. Critical Behaviour in Finite Temperature QCD 

We want to consider here in particular the role of the 'bare' quark mass m q in the 
Lagrangian 

C = -\f; v F^ v - ^(i^-g^A^p -m q ^ a (1) 
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for the critical behaviour of statistical QCD. From the point of view of pure gauge theory, 
the introduction of m q breaks explicitly the center symmetry of the Lagrangian and 
thus is like turning on an external magnetic field H ~ l/m q in the corresponding 
spin theory; for infinite quark mass, we recover the thermodynamics of pure gauge theory, 
just as H = leads back to the Z^-symmetric spin system. For sufficiently small quark 
masses, one might thus expect a smooth temperature variation of the Polyakov loop ||, 
similar to the magnetization pattern of spins in a strong external field. However, such 
a behaviour is not observed; even for quark masses approaching the chiral limit, the 
Polyakov loop 

L(T,m q )~eM-V(m q )/T} (2) 

varies very rapidly at those 'critical' temperatures at which energy density and chiral 
condensate vary strongly; here V(m q ) denotes the interquark potential in the limit of large 
separation. The reason for this persistence of an 'almost' critical behaviour of L(T, m q ) in 
T seems evident: the spontaneous breaking of the chiral symmetry for m q = produces 
an effective 'constituent' quark mass M q ~ Mh/2, equal to about half the mass Mh of 
the basic (non-Goldstone) meson. For m q ^ 0, chiral symmetry is explicitly broken and 
M q {m q ) thus varies with m q in the form shown in Fig. 2. Equivalently, the potential string 
breaks in the large distance limit not for 2m q , but rather for ~ 2M q , so that even for 
m q — > 0, a finite energy is still required for string breaking. It is M" 1 , not the inverse 
bare quark mass m" 1 , which acts as external field, and hence this field always remains 
quite weak. The range of variation of L(T,m q ) thus lies between the singular pure gauge 
form L(T, m q = oo) and a limiting form L(T, m q = 0) which still varies strongly with T 
(Fig. 3). ' 




Figure 2: Variation of the constituent quark mass with bare quark mass. 



The thermal pattern of L(T, m q ) should be contrasted to that of the chiral condensate 
(^■^(T, m g ). For m q = 0, the spontaneous breaking of the chiral symmetry at low 
temperatures disappears when T = T x , and from there on chiral symmetry is restored. 
The introduction of a non- vanishing bare quark mass m q explicitly breaks chiral symmetry, 
so that here m q (and not its inverse) acts like an external field in the corresponding 
0(4) spin system ||; only for m q = is there genuine spontaneous symmetry breaking, 
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Figure 3: The temperature dependence of the Polyakov loop (a) for pure gauge theory 
(m q — > oo), (b) for QCD with light quarks (0 < m q < oo), and (c) in the chiral limit 

K = o). 

and when m q becomes large, the temperature variation of {ipip)(T,m q ) is washed out 
completely. 

The 'transition' line between m q = oo and m q = in Fig. 1 can thus be character- 
ized as follows: the singular behaviour of L(m q = oo) changes into an 'almost' singular 
behaviour as the quark mass becomes finite, with L(M q [m q = 0]) as the chiral limit. Any 
(spontaneous or explicit) breaking of chiral symmetry assures that M q ^ and hence lim- 
its the strength of the effective external field. The associated 'almost' singular behaviour 
of thermal observables such as entropy or energy density is thus related to deconfinement 
as far as the overall variation is concerned; the chiral features determine how 'smooth' 
this variation can become, and explicit or spontaneous chiral symmetry breaking keeps it 
always very pronounced. - While these considerations provide a reason why the decon- 
finement transition, as defined by the behaviour of the Polyakov loop, 'almost' persists in 
the chiral limit, they do not explain why deconfinement and chiral symmetry restoration 
effectively coincide at zero baryon number density; all arguments appear to remain valid 
if chiral symmetry were restored above deconfinement. 

We now turn to percolation theory, with the aim of obtaining a characterization of 
deconfinement which remains valid for all values of the quark mass m q . 

3. Percolation and the Ising Model 

For simplicity, consider a two-dimensional square lattice of linear size L; we randomly 
place identical objects on iV of the L 2 lattice sites. With increasing N, adjacent occupied 
sites will begin to form growing clusters or islands the sea of empty sites. Define n p to 
be the lowest value of the density n = N/L 2 for which the origin belongs on the average 
(i.e., calculated by averaging over many randomly generated configurations) to a cluster 
reaching the edge of the lattice. In the limit of infinite lattice size, we then have 

(rr \ ftp 
, n>n p , (3) 
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for n approaching n p from above, where P[n) denotes the probability that the origin 
belongs to an infinite cluster. Since P(n) = for all n <n p and non-zero for all n > n p , 
it constitutes an order parameter for percolation: f3 p = 5/36 is the critical exponent 
which governs the vanishing of P(n) at n = n p in two dimensions; in three dimensions, it 
becomes (3 = 0.41 §. 

For n < n p , a quantity of particular interest is the 'mean' cluster size S(n), defined 
as the average size (the number of connected occupied sites) of a cluster containing the 
origin of the lattice. As n approaches n p from below, S(n) diverges at the percolation 
point as 

s d (n) ~ (n p - n)~ 7p , n < n p (4) 

with 7 p = 43/18 (1.80) as the d — 2 (3) critical exponent for the divergence. 

Instead of the site percolation just described, a random placement of bonds on the 
lines between adjacent sites leads to bond percolation, with the same critical exponents. 
Clusters now consist of regions of connected bond lines. While the critical behaviour 
is universal, the value of the percolation threshold n p depends on the type of percola- 
tion considered; on a two-dimensional square lattice, it is 0.5 for bond and 0.59 for site 
percolation ||. 

We now turn to the Ising model, defined by the Hamiltonian 

Hi = -Jj2 s i s j ~ H ^Z s h Si = ±1, (5) 

where the first sum runs only over nearest neighours on the lattice, with J denoting the 
exchange energy between spins and H an external field. For H = 0, the Hamiltonian (3) 
has a global Z2 invariance (sj — > —Si V i), and the magnetization m = (s) probes whether 
this invariance is spontaneously broken. As is well-known, such spontaneous symmetry 
breaking occurs below the Curie point T c , with 

m (T,H = 0)~(l-^ (6) 

governing the vanishing of m(T, H = 0) as T — > T c from below. The well-known Onsager 
solution gives with (3 m = 0.125 for d = 2 a 10 % smaller value than the (3 P = 5/36 ~ 0.139 
found for the percolation exponent. 

Since the Ising model also produces clusters on the lattice, consisting of connected 
regions of aligned up or down spins, the relation between its thermal critical behaviour at 
T c and the onset of geometric percolation is an obvious question which has been studied 
extensively in recent years. Two important distinctions have emerged: the nature of what 
is defined as a cluster J7[, and the behaviour of clusters in a non- vanishing external field 
1,1- 

The geometric clusters considered in percolation theory consist simply of connected 
regions of spins pointing in the same direction. In the Ising model, there is a thermal 
correlation between spins on different sites; this vanishes for T — ► 00. Correlated regions 
in the Ising model (we follow the usual notation and call them 'droplets', to distinguish 
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them from geometric clusters) thus disappear in the high temperature limit. In contrast, 
the geometric clusters never vanish, since the probability for a finite number of adjacent 
aligned spins always remains finite; it increases with dimension because the number of 
neighbours does. Hence from the point of view of percolation, there are more and bigger 
clusters than there are Ising droplets. As an immediate consequence, percolation occurs 
in general before the onset of spontaneous magnetization. 

If percolation is to provide the given thermal critical behaviour, the definition of cluster 
has to be changed such that the modified percolation clusters coincide with the correlated 
Ising droplets @. This is achieved by assigning to some pairs of adjacent aligned spins in 
a geometric cluster an additional bond correlation, randomly distributed with the density 

n b = l- 6' 2J / fcT , (7) 

where [—2J/kT] just corresponds to the energy required for flipping an aligned into an 
opposing spin. The modified percolation clusters then consist of aligned spins which are 
bond connected. Only for T = are all aligned spins bonded; for T > 0, some aligned 
spins in a purely geometric cluster are not bonded and hence do not belong to the modified 
cluster or droplet. This effectively reduces the size of a given geometric cluster or even 
cuts it into several modified clusters. For T — > oo, rib — * 0, so that the geometric clusters 
still in existence there are not counted as droplets, solving the problem mentioned in the 
previous paragraph. 

With such a superposition of site and bond percolation ('s/b'), full agreement between 
percolation and thermal critical behaviour of the Ising model is achieved. An infinite 
percolation cluster is now formed for the first time at T c , the cluster size coincides with 
that of the correlated regions in the Ising model, and numerical simulations show that 
the effective critical exponents for the new P s /b{n) or S s /b(n) become those of the Ising 
model. 

However, for H ^ 0, Ising model and percolation theory results do not coincide, nei- 
ther for purely geometric nor for site/bond defined percolation clusters. The Ising model 
for H ^ does not lead to any singularity as function of T and hence does not show any 
critical behaviour ]T(| ; the symmetry responsible for the onset of spontaneous magne- 



tization is now always broken and m(T, H ^ 0) ^ for all T. On the other hand, the size 
of connected regions of aligned and bonded spins increases with decreasing temperature, 
and above some critical temperature it diverges. Hence percolation will occur for any 
value of if. At H — oo, all spins are aligned, leaving the bonds as the relevant variables; 
the system thus percolates at the critical density for bond percolation, which is 1/2 in two 
and 0.249 in three dimensions, for square and cubic lattices, respectively. From Eq. (^|), 
this corresponds to temperatures kT^/J = 2.89 and 6.99 in two and three dimensions, 
to be compared to the Curie temperatures kT c / J = 2.27 and 4.51, respectively. The 
corresponding values at finite H lie between T c at H = and T& at H = oo and define the 
so-called Kertesz line || ||; see Fig. 4. What happens at this line in terms of dynamics? 

Since the Ising model does not lead to singular behaviour in T for H ^ 0, in this case 
even the modified percolation picture does not agree with the Ising result. The reason 
for this is physically quite evident. At sufficiently high temperatures, only finite size 



6 




Figure 4: The Kertesz line, defining the percolation temperature in the Ising model as 
function of the external field H; T c denotes the Curie temperature, T b the pure bond 
percolation temperature. 



s/b clusters exist, so that the percolation probability vanishes; at low temperatures, on 
the other hand, the s/b cluster size diverges, and hence P(n) remains as genuine order 
parameter. The existence of a non- vanishing external field, however, partially aligns in its 
direction the spins of the different disconnected finite size clusters at high temperature U, 
so that the overall magnetization never vanishes. As a result, the critical behaviour due to 
percolation persists, while that related to the thermal magnetization pattern disappears. 

As is well-known, the Ising model can be reformulated as a lattice gas, where its critical 
behaviour describes in a rudimentary way the vapour-liquid transition. The non-singular 
behaviour for H ^ then translates into a continuous cross-over from gas to liquid for 
sufficiently high pressure; the two phases are assumed to become indistinguishable in 
the relevant pressure-temperature regime. However, the percolation effect discussed in 
the previous paragraph seems to contradict the established credo that 'nothing happens 



above the critical point'. It is observed H, ITTJ that for temperatures below the Kertesz 



line, Ising droplets have a finite surface tension; this vanishes along the Kertesz line, 
and on its high temperature side, the droplet energy becomes proportional to their bulk 
energy. Note that such a behaviour does not contradict the non-singular thermodynamic 
behaviour of the model, since the variable that changes along the Kertesz line, the surface 
tension of droplets, is not expressible in terms of thermodynamic observables. 

In the same vein, we note another illustration of the relation between percolation 
and thermal critical behaviour. Consider bond percolation in the two-dimensional Ising 
model and imagine that current can flow between two or more bonded sites. In this 
case, conductivity sets in at the percolation point, independent of the thermal critical 
behaviour of the Ising model; the system is non-conducting below the percolation point 
and conducting above it. This remains true also in s/b percolation and is independent of 
whether clusters are partially aligned by an external field. In other words, we have two 
independent critical phenomena, the onset of conductivity and the critical behaviour of 
Ising thermodynamics. This suggests a possible framework for a general characterization 
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of deconfinement. Before addressing this, however, we want to point out some aspects 
encountered in the extension of percolation to continuum systems @. 

4. Continuum Percolation 

We start again in two dimensions and consider the distribution of discs of radius r 
over some plane area, allowing overlap; as example, we take a circle of radius R ^> r 
(see Fig. 5). The average cluster size s c \ and density n c \ (number of discs per cluster) 
is determined also here by averaging clusters containing the origin over many randomly 
generated configurations. These quantities are studied as function of the overall density n, 
i.e., the number of discs per circle ttR 2 . In the limit R —>■ oo, the percolation probability 
P{n) and the divergence of the cluster size S{n) are then given by Eqs. @) and @), with the 
same critical exponents; the same remains true for the corresponding three-dimensional 
case, with overlapping spheres instead of discs. 




Figure 5: Two-dimensional continuum clustering. 

In the continuum, the percolation threshold becomes || |T^j 

1.18 



7T7* 2 

for d = 2 and 

n ~ ° M (9) 

for d = 3. In both cases, there is considerable overlap, and at the percolation point, the 
fraction of the area (volume) covered by discs (spheres) is 1 — exp(— r] p ), where r] p ~ 1.18 
(0.34) for d = 2 (3). 

Besides the discs and spheres considered here, the percolation of cylinders or ellipsoids 
will turn out to be of interest to us. While the critical behaviour at the percolation 



8 



threshold is universal, i.e., dependent only on the space dimension of the problem, the 
threshold value itself is affected by the geometry of the objects distributed in space. 
For ellipsoids, percolation has been studied numerically in considerable detail [I3|: as 
representative illustration, the threshold density for percolation is found to be 

"> " WWh^ (10) 

in the case of an aspect ratio a/b = 4 (a is the long, b the short axis of a prolate ellipsoid). 

In section 3, we had noted the Coniglio-Klein cluster definition 0, which resulted in 
the coincidence of percolation and thermal critical behaviour of the Ising model. For our 
later considerations, the extension of this definition to the continuum is necessary. This, 
however, for the moment appears to remain an interesting open problem in percolation 
theory fll4] . 



5. Percolation in QCD 

We now turn to the application of percolation in the study of strongly interacting 
matter and consider first pure colour SU (N) gauge theory. Here the problem is straight- 
forward: we have to extend the Coniglio-Klein approach to continuum percolation, 
i. e., construct a cluster definition which moves the critical behaviour at the percolation 
threshold from the universality class of geometric percolation to that of thermal critical 
behaviour of a Zn spin system ]T5| . If this is possible, glueball percolation will lead to 



the correct deconfinement behaviour as obtained in terms of the Polyakov loop as order 
parameter. 

For the case of QCD with dynamical quarks, it had already been noted in the first 
studies Jl], HJ that the percolation of spheres of hadronic size cannot mean deconfinement: 
using Eq. (|]) and a hadron radius r = 1.0 fm, the percolation threshold becomes n p ~ 0.08 
fm -3 , i.e., it lies below normal nuclear density. Hence at this density, the quarks must 
still remain effectively localized in colour singlet states, so that the percolation point here 
corresponds to the formation of connected nuclear matter. Deconfinement, it was argued, 
occurs only when percolation sets in for some smaller 'hard core' of hadrons. The question 
was how to define this. 

In QCD, two fundamental triplet charges are connected by a colour flux tube or string, 
with a ~ 0.16 GeV 2 as string tension. The transverse radius r of these strings has the 



form [jlB] , PH 



r 2 = — ln(/// ), (11) 

7T<7 



where Iq ~ 0.1 — 0.3 fm is the string formation length [IR |T9 . Extensive lattice QCD 
studies find that in the range 0.5 < / < 2 fm, r ~ 0.2 — 0.3 fm and varies at most weakly 
TT|, [HI. As a result, a hadron appears as a cylindrical string of diameter 2r and length 



/ ~ 1 fm, oscillating about its center of mass to fill a sphere of that radius. While the 
energy density in the string of length I and radius r is of order 5 GeV/fm 3 , the oscillation 
spreads this over a sphere of radius I and thus reduces the energy density to the 0.2 - 0.3 
GeV/fm 3 observed for a normal (non-Goldstone) meson. 
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Given such a partonic substructure of a hadron, it becomes clear that deconfinement 
requires the percolation of the underlying strings. Evidently (compare Eqs. (§) and (|1C|) ) 
this leads to a much higher percolation threshold: while spheres of radius r = 1 fm 
percolate at ^sphere ^ g gg f m -3^ threshold for strings of that length and of radius 
0.25 fm (taken as corresponding ellipsoids) becomes more than ten times higher, with 
n string ^ ^ gg f m -3 p or (non-Goldstone) hadrons having a generic mass of about 1 GeV, 
this leads to deconfinement at an energy density around 1 GeV/fm 3 , in general agreement 
with lattice results. Equating Eqs. @ and ( |10"1) gives 

r core ~ 0.46 fm (12) 

as the radius for the percolation of spheres equivalent to the given string percolation. 
The 'hard core' needed to define deconfinement thus has a radius of about half that of a 
hadron. 

Once the problem of extending the Coniglio-Klein cluster definition to the continuum 
is solved, we can propose to identify the onset of deconfinement with the percolation 
threshold for strings of radius r and hadronic length I. In the limit of infinite quark mass, 
the critical behaviour of SU(N) gauge theory and the corresponding percolation then 
coincide by virtue of the Coniglio-Klein construction. For any finite value of the bare 
quark mass, the QCD counterpart of the Kertesz line will separate a high density colour 
conducting deconfined region - the percolation regime - from a low density confined colour 
insulator, the finite cluster regime. In the confined region, the average Polyakov loop does 
not vanish; it is only exponentially small, L ~ exp(— M^/T), because of occasional thermal 
string breaking. On the other hand, the average cluster size in this region is finite, so 
that the percolation probability remains zero up to the percolation threshold. 

6. Deconfinement in Nuclear Collisions 

While deconfinement in an equilibrated QCD medium thus should start at the perco- 
lation threshold for a system of randomly oriented strings, the pre-equilibrium situtation 
encountered in nuclear collisions provides an alignment or polarization of these strings. 
Consider a central high energy collision of two identical heavy nuclei, in which essentially 
all nucleons undergo several interactions. These occur in such rapid succession that the 
nucleons cannot 'recover' before being hit again. Each individual nucleon-nucleon colli- 
sion establishes a colour flux tube or string between the collision partners; the additive 



quark model [^0| suggests that at present (SPS) energies it connects two triplet colour 
charges, so that it this just the string considered above [T5|, O. For nuclear collisions 
in their early stage, we thus obtain a spaghetti-like structure of intertwined, overlapping 
QCD strings, and a cut in the transverse plane results in the picture shown in Fig. 5, 
giving a distribution of transverse string areas over a circle of nuclear area. Hence the 
two-dimensional continuum percolation theory presented above is the tool to study what 
happens when these strings overlap more and more to form connected spatial regions of 
increasing size in the transverse planeQ. One important modification is that the overall 
size of the transverse area is now fixed by nuclear dimensions, requiring the study of finite 
size cluster formation. 



2 The first study of this kind was carried out in a specific phenomenological string model plj 
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We first consider the cluster density in the transverse plane, then the geometric cluster 
size as function of the density of the overall density of string discs. To have a specific 
case, we set R/r = 20, corresponding to R — 5 fm and r = 0.25 fm; to eliminate 
some of the dependence on the specific choice of parameters, it is convenient to make all 
quantities (n, n c \ and s c \) dimensionless by defining them in terms of irr 2 . In Fig. 6 the 
resulting cluster density n c \ is seen to increase monotonically with the overall density n. 
As expected from percolation theory, the size s c \ of the cluster, in contrast to its density, 
shows a dramatic variation with n. As noted in section 4, s c \ diverges for R/r —>■ oo at the 
percolation point n p = 1.175, i.e., when the overall density is somewhat above one disc 
per disc area. Of course there is a considerable overlap of discs, and at n p , the ratio of 
the area covered by discs to the overall area becomes 1 — exp[— 1.12] ~ 0.67; it approaches 
unity only for n — > oo. For finite R/r, the critical behaviour is softened by finite size 
effects, and for the value R/r = 20 chosen above, we obtain the cluster size behaviour 
shown in Fig. 7. It indeed shows the strongest variation around the percolation point. 
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Figure 6: Cluster density (a) and cluster size (b) as function of overall density, in two 
dimensions, with r = 0.25 fm, R = 5 fm; n p denotes the infinite area percolation point. 



Deconfinement is expected to set in when there is enough 'internetting' between inter- 
acting nucleons P3, i.e., when the collision density within a cluster becomes sufficiently 
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high; we denote this value by n c cl . It is reached at a certain value n c = n{n c c ^) of the overall 
density, which can, but does not need to be the ultimate percolation point. By Fig. 7, 
this identifies a corresponding critical cluster size s c cl = s c i(n£i). In other words, requiring 
a critical collision density automatically forces the cluster to have a certain minimum size. 
This result answers the question of how much medium is needed before one can speak 
about a new phase: systems of geometric size s c \ < will on the average not have a 
string density sufficient for deconfinement. 

The onset of deconfinement in nuclear collisions is thus governed by two parameters: 
the deconfinement string density n c cX (which could be the percolation value) and the trans- 
verse string size r (which is roughly determined by QCD |Tj| [17j). The size of the nuclear 
transverse area irR 2 is determined by nuclear geometry. With n c d and r fixed, we can then 
determine the size s c cl of the bubbles of deconfined medium present at the deconfinement 



threshold. Realistic nuclear distributions 22 and a Glauber-based formalism must be 



used to calculate the distribution of collisions in the transverse plane as function of the 



impact parameter in A — B collisions [23]. The distribution of the strings in the transverse 
plane is then no longer random, but governed by the nucleon-nucleon collision density thus 
determined. The resulting formalism can be used to study the onset of deconfinement, 
using the predicted J/ip suppression as indicative signal J24J. Here the finite cluster size 



at the deconfinement threshold leads to an abrupt onset of such suppression |25| , which 
agrees well with the recently reported anomalous J / ip suppression in Pb — Pb collisions 
at the CERN-SPS pfj. Moreover, the critical density n c needed here also agrees rather 
well with the (infinite volume) percolation point. 
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